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SIGNAL  PROCESSING  FAULT  DETECTION  SYSTEM 


STATEMENT  OF  GOVERNMENT  INTEREST 
[0001]  The  invention  described  herein  may  be  manufactured  and 

used  by  or  for  the  Government  of  the  United  States  of  America 
for  governmental  purposes  without  the  payment  of  any  royalties 
thereon  or  therefore. 

CROSS  REFERENCE  TO  OTHER  PATENT  APPLICATIONS 
[0002]  None. 

BACKGROUND  OF  THE  INVENTION 
(1)  FIELD  OF  THE  INVENTION 

[0003]  The  present  invention  is  directed  to  the  detection  and 
analysis  of  cracks,  asymmetries,  and  imperfections  within  a 
material  structure.  In  particular,  the  present  invention  is 
directed  to  a  diagnostic  tool  that  determines  the  linear  and 
nonlinear  characteristics  of  measured  vibration  signals  in  order 
to  use  the  most  appropriate  signal  processing  technique  for  the 
identification  of  faults  in  structures  and  to  determine 
quantitatively  the  predisposition  of  a  structure  to  material 
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(2)  DESCRIPTION  OF  THE  PRIOR  ART 

[0004]  Dynamic  loading,  such  as  vibration,  can  be  recorded  by- 
strain  sensors.  The  strain  sensor  output  signals  in  the  form  of 
time  series  data  can  be  analyzed  to  determine  the  condition  of  a 
physical  structure.  Sensor  signals  from  a  "healthy"  structure 
(i.e.  a  structure  free  of  cracks,  asymmetries  and  imperfections) 
are  typically  linear.  The  sensor  signals  from  a  structure  with 
fatigue  cracks  and  or  other  imperfections  are  typically 
nonlinear.  The  nonlinearity  of  time  series  data  from  sensors 
recording  vibrations  in  structures  with  fatigue  cracks  can  be 
classified  as  a  discontinuity  or  as  the  simultaneous  emission  of 
several  frequencies.  With  the  advent  of  increased  computational 
power  in  computers,  there  exists  the  computing  capability  to 
perform  the  calculations  necessary  to  apply  non-linear 
mathematical  analysis  techniques  to  time  series  data  obtained 
from  sensors,  such  as  strain  sensors. 

[0005]  There  are  several  mathematical  analysis  techniques  to 
determine  the  linear  and/or  nonlinear  characteristics  of  strain 
sensor  signals  which  represent  the  vibrations  in  a  structure. 

One  such  technique  is  mutual  information  analysis.  Mutual 
information  analysis  is  a  mathematical  analysis  technique 
derived  from  Information  Theory.  When  dealing  with  random 
variables,  a  mutual  information  analysis  will  seek  to  determine 
the  amount  of  information  that  one  random  variable  contains 
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about  another  random  variable  and  vice  versa.  This  sort  of 
determination  serves  as  a  measure  of  dependence  between  the 
first  and  second  random  variables.  Mutual  information  analysis 
can  be  used  to  test  the  dependencies  between  two  sets  of  time 
series  data,  such  as  the  data  obtained  from  strain  sensors. 

[00 OS]  Another  mathematical  analysis  technique  to  determine 

the  linear  and  nonlinear  characteristics  of  strain  sensor 
signals  is  wavelet  analysis  which  is  a  linear  mathematical 
analysis  technique  that  can  analyze  discontinuities  and  edge 
effects.  Wavelet  analysis  can  define  wavelets  in  either  the 
real  domain  (referred  to  as  real  wavelets)  or  the  complex  domain 
{referred  to  as  analytic  wavelets) .  Real  wavelets  are  suitable 
for  identifying  discontinuities  and  data  compression.  Analytic 
wavelets  are  suitable  for  capturing  frequency  content  within  a 
signal  and  therefore  isolating  simultaneous  frequency  emissions. 
[0007]  Other  mathematical  analysis  techniques  to  determine 
the  linear  and  nonlinear  characteristics  of  strain  sensor 
signals  in  the  form  of  time  series  data  are  the  surrogate  data 
method,  fast  Fourier  transforms  and  the  phase  space  method. 

These  methods  are  utilized  (sometimes  with  other  information  or 
signal  processing  diagnostic  techniques)  to  ascertain  the  linear 
or  nonlinear  characteristics  of  the  measured  data. 

[0008]  The  above-described  mathematical  analysis  techniques 
are  capable  of  detecting  the  linear  and/or  nonlinear  "behavior" 
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of  cracks,  asymmetries,  and  imperfections  within  a  material 
structure.  In  the  past,  methods  for  evaluating  cracks  and 
imperfections  in  material  structures  have  relied  primarily  on 
linear  mathematical  analysis  techniques  as  applied  to  time 
series  data  obtained  from  strain  sensors  (i.e.  linear  signal 
processing  techniques) ,  rather  than  non-linear  mathematical 
analysis  techniques  (i.e.,  non-linear  signal  processing).  Some 
linear  signal  processing  techniques,  however,  fall  short  of 
identifying  material  failures  or  incongruities  where  such 
failures  or  incongruities  are  due  to  "nonlinearity"  in  the 
failure  generating  mechanism (s) .  In  order  to  truly  capture  the 
range  of  potential  material  failures  and  discontinuities  in  a 
material  structure,  nonlinear  signal  processing  techniques  must 
also  be  considered  to  specifically  determine  any  nonlinear 
behavior  in  the  material  structures.  Determining  the  proper 
approach  and  combination  of  linear  and  nonlinear  data  signal 
processing  techniques  for  evaluating  structures  in  an  efficient 
manner  has  become  the  challenge. 

SUMMARY  OP  THE  INVENTION 

[0009]  It  is  a  general  purpose  and  object  of  the  present 
invention  to  efficiently  obtain  linear  and  nonlinear 
characteristics  from  time  series  data  obtained  from  a  sensor 
system  in  order  to  analyze  the  linear  and  nonlinear  behavior  of 
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cracks,  asymmetries,  and  imperfections  within  a  material 
structure,  and  thereby  detect  such  flaws  in  the  structure. 

[0010]  It  is  a  further  object  of  the  present  invention  to 
determine  the  proper  approach  and  combination  of  linear  and 
nonlinear  signal  processing  techniques  for  efficiently 
evaluating  time  series  data  obtained  from  strain  sensors 
attached  to  material  structures. 

[0011]  The  above  objects  are  accomplished  with  the  present 
invention  by  a  system  employing  a  series  of  strain  sensors  and  a 
data  processor  that  implements  a  method  to  preliminarily  analyze 
the  linear  and  nonlinear  characteristics  of  the  time  series  data 
generated  by  the  strain  sensors  as  they  sense  vibrations  in  the 
structure,  and  then  determine  whether  further  linear  or 
nonlinear  analysis  on  the  time  series  data  is  warranted.  After 
further  analysis  is  performed  on  the  time  series  data,  the 
resulting  values  are  compared  to  baseline  values  to  determine  if 
there  are  any  significant  imperfections  in  the  material 
structure  that  would  indicate  a  damaged  structure.  A  graphical 
user  interface  joined  to  the  data  processor  is  used  to  alert  a 
user  of  potential  damage  in  the  material  structure . 


BRIEF  DESCRIPTION  OF  THE  DRAWINGS 
[0012]  A  more  complete  understanding  of  the  invention  and 
many  of  the  attendant  advantages  thereto  will  be  more  readily 
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appreciated  by  referring  to  the  following  detailed  description 
when  considered  in  conjunction  with  the  accompanying  drawings, 
wherein  like  reference  numerals  refer  to  like  parts  and  wherein: 
[0013]  FIG.  1  is  a  block  diagram  of  the  system  components 
including  the  sensors,  the  data  processor,  and  the  user 
interface ; 

[0014]  FIG.  2  is  a  flow  chart  of  the  method  used  in  the  fault 
detection  system; 

[0015]  FIG.  3  is  an  illustration  of  a  helicopter  flexibeam 
used  as  a  test  structure; 

[0016]  FIG.  4  is  a  graph  of  the  time  series  data  representing 
the  output  of  the  strain  sensors  attached  to  both  appendages  of 
the  flexibeam; 

[0017]  FIG.  5  is  a  graph  of  the  mutual  information  analysis 
of  the  "healthy"  appendage; 

[0018]  FIG.  6  is  a  graph  of  the  mutual  information  analysis 
of  the  damaged  appendage; 

[0019]  FIG.  7  is  a  three  dimensional  graph  of  the  wavelet 
transmittance  function  of  both  appendages; 

[0020]  FIG.  8  is  an  illustration  of  a  damaged  and  a  "healthy" 
cantilever  beam  both  having  sensors  attached  to  each;  and 
[0021]  FIG.  9  is  a  graph  of  the  mutual  information  analysis 
of  the  sensor  data  obtained  from  both  cantilever  beams. 
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DETAILED  DESCRIPTION  OF  THE  INVENTION 
[0022]  Referring  initially  to  FIG.  1,  there  is  shown  a  block 
diagram  of  the  signal  processing  fault  detection  system  10  of 
the  present  invention  consisting  of  an  array  of  sensors  15 
capable  of  detecting  vibrations  in  a  structure,  a  data  processor 
20,  and  a  graphical  user  interface  25.  Sensors  15  such  as 
piezoelectric  strain  sensors  are  deployed  along  a  physical 
structure  with  the  ultimate  goal  of  assessing  the  condition  of 
the  structure.  Other  sensors  that  are  capable  of  detecting 
vibrations  could  also  be  used.  The  sensors  15  are  connected  to 
the  data  processor  20,  such  as  a  computer  system,  that  can 
receive  sensor  data  and  perform  linear  or  nonlinear  signal 
processing  on  the  sensor  signals. 

[0023]  Referring  to  FIG.  2  there  is  shown  a  flow  chart  of  the 
method  to  determine  linear  and  nonlinear  characteristics  of  time 
series  data  30  obtained  from  the  array  of  strain  sensors  15 
deployed  on  a  physical  structure  to  detect  and  record  vibrations 
in  the  structure.  The  first  step  40a  in  the  flow  chart  is  to 
create  a  computer  model  of  the  physical  structure  that  is 
relatively  free  of  structural  faults.  Such  a  model  can  be 
created  using  standard  modeling  software  well  known  in  the  art. 
The  computer  model  will  provide  time  series  data  representing 
vibration  signals  that  will  serve  as  a  baseline  for  the 
structure.  Alternative,  the  first  step  40b  is  to  obtain  a 
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baseline  of  time  series  data  from  strain  sensors  15  disposed 
about  an  actual  physical  structure  that  is  relatively  free  of 
structural  faults. 

[0024]  The  second  step  in  the  flow  chart  45a  is  to  quantify 
the  vibration  signals  obtained  from  the  computer  modeled 
structure  and  then  perform  a  mutual  information  analysis  on  the 
signals  to  obtain  a  baseline  of  nonlinear  values. 

Alternatively,  the  second  step  45b  is  to  quantify  the  actual 
historical  data  from  the  physical  structure  and  then  perform  a 
mutual  information  analysis  on  the  signals  to  obtain  a  baseline 
of  nonlinear  values.  The  baseline  serves  the  purpose  of 
establishing  a  threshold  for  comparison.  For  example,  when 
measured  quantified  vibration  signals  are  obtained  from  sensors 
disposed  over  similar  structures  or  the  same  structure,  and  the 
quantified  vibration  signals  are  compared  to  and  exceed  the 
baseline,  then  it  is  an  indication  that  there  are  faults  in  the 
similar  structures  or  the  same  structure. 

[0025]  The  third  step  in  the  flow  chart  50  is  to  obtain 
experimental  time  series  data  of  vibration  signals  from  an  array 
of  strain  sensors  15  deployed  on  a  similar  structure  or  the  same 
structure  for  the  purpose  of  detecting  any  material  faults  in 
the  similar  structure  or  the  same  structure. 

[0026]  The  fourth  step  in  the  flow  chart  55  requires  a 
preliminary  mutual  information  analysis  be  performed  on  the 


8 


Attorney  Docket  No,  97588 


experimental  time  series  data  to  determine  whether  the  data  has 
linear  or  nonlinear  characteristics.  The  degree  that  data  has 
linear  or  nonlinear  characteristics  requires  a  general  and  not 
specific  determination  because  in  theory  each  series  of  data  can 
represent  a  unique  physical  situation.  There  can  be  a  fault  in 
a  structure  that  exhibits  nonlinear  behavior  but  the  degree  of 
nonlinear  behavior  does  not  exclude  the  possibility  of  using  a 
linear  analysis  technique  to  detect  the  fault. 

[0027]  Equation  1  is  the  general  mutual  information  equation 
used  in  the  fourth  step  55  of  the  method  illustrated  in  FIG  2 . 
The  general  mutual  information  equation  is  a  probabilistic 
equation  used  here  to  assess  the  information  between  two  random 
variables . 


rA.a  = 


PamW 

L  PM)P.(*>j)} 


(1) 


[0028]  PAa(atbj)  is  the  joint  probability  of  a  set  A  =  {at}  and  a 
set  B  =  y>l\,  while  (a,)  and  PB(bj)  are  the  marginal  individual 

probabilities  of  each  set  of  A  and  B.  For  example,  set  B  can 
be  a  time-delayed  version  of  set  A  . 

[0029]  The  fifth  step  in  the  flow  chart  60  is  a  decision  step 
that  requires  a  determination  of  whether  the  time  series  data 
has  nonlinear  characteristics  based  on  the  dependencies  between 
the  time  series  data  sets  A  and  B.  If  set  A  and  set  B  have 
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dependencies,  then  that  is  an  indication  that  the  time  series 
data  has  linear  characteristics  because  there  is  a  correlation 
between  the  signals.  If  the  time  series  data  has  linear 
characteristics,  then  a  linear  signal  processing  analysis  will 
be  performed.  If  the  time  series  data  does  have  nonlinear 
characteristics,  then  further  nonlinear  signal  processing 
analysis  will  be  performed. 

[0030]  According  to  the  sixth  step  in  flow  chart  65a,  when 
the  data  is  determined  to  have  linear  characteristics  as 
described  above,  a  linear  analysis  such  as  wavelet  transform 
analysis  as  expressed  in  Equation  2  or  fast  Fourier  transform 
(FFT)  analysis  as  expressed  in  Equation  3  is  performed.  The 
wavelet  transformation  equation  is  described  according  to: 


(2) 


where  f(t)  is  the  measured  function  (i.e.  the  time  series  data 
captured  from  the  sensors)  and  y/"»,u  is  identified  as  the  complex 

conjugate  of  the  mother  wavelet.  The  variable  1  represents  time, 
a  is  defined  as  the  scaling  parameter,  and  u  is  termed  the 
translation  within  the  time  domain. 

[0031]  The  Discrete  Fast  Fourier  Transform  is  written  as 


(3) 
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where  x(t„)  is  the  input  signal  from  the  sensors  at  time  t„  usually 
in  (sec)  ,  t„  is  the  nth  sampling  instant  in  (sec)  .  X((OfJ  is  the 
spectrum  of  a-  at  frequency  ( 0 *  ,  ftf*  is  the  kth  frequency  sample  in 

(rad  /sec) ,  and  N  is  the  number  of  samples  and  number  of 
frequency  samples. 

[0032]  According  to  the  seventh  step  of  the  flow  chart  70a, 
the  results  of  the  linear  analysis  of  the  time  series  data  are 
compared  to  the  baseline  values  obtained  in  the'  second  step  45a 
or  45b  to  approximate  the  damage  level,  if  any,  to  the 
structure.  The  level  of  damage  is  an  approximation  based  on  the 
differences  between  the  analyzed  time  series  data  and  the 
baseline  values.  The  magnitude  of  the  difference  is  directly 
related  to  acceptable  tolerances,  which  are  based  on  the 
application  of  the  structure  itself.  A  static  structure  may 
allow  for  more  of  a  difference  between  the  analyzed  data  and 
baseline  values  than  a  dynamic  structure. 

[0033]  According  to  the  eighth  step  of  the  flow  chart  75a  a 
decision  must  be  made  as  to  whether  to  generate  an  alert  because 
the  structure  is  seriously  damaged.  If  there  is  a  serious  level 
of  damage  as  indicated  by  a  large  difference  beyond  the 
acceptable  baselines,  then  the  graphical  user  interface  will 
create  a  visual  alert  80a  for  a  user.  If  there  is  not  a  serious 
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level  of  damage,  the  next  set  of  time  series  data  of  vibration 
signals  is  analyzed  in  the  final  step  85a. 

[0034]  Alternatively,  at  the  sixth  step  of  the  flow  chart 
65b,  after  the  data  is  determined  to  be  nonlinear  as  described 
above,  a  further  nonlinear  analysis  is  performed  on  the  time 
series  data.  Nonlinear  analysis  techniques  include  but  are  not 
limited  to  mutual  information  as  described  in  Equation  1, 
combined  with  surrogate  data  analysis  as  illustrated  in  Equation 
4  or  phase  space  analysis  as  is  also  known  in  the  art. 

[0035]  The  surrogate  data  method  where  the  phase  is 
randomized  is  expressed  in  Equation  4: 

(  J(")exP['V(»)]  for  n  =  l:Af/2-l  ] 

[s(«)exp[-(^(W-rt),/br  n  =  Nil  +  1 :  N -1 J 


[003  6]  Sequence  {Sfal}  is  the  overall  phase -randomized  FFT 
signal,  and  sequence  { <j>{n) } ,  n=\:N/l-\,  is  randomly  shuffled,  over 
2n  radians,  from  the  original  time  series,  is  the  complex 

amplitude  of  the  FFT  of  the  original  time  series.  In  the 
surrogate  method,  the  inverse  FFT  is  calculated  for  the  modified 
and  this  transformed  time  series  is  called  the  surrogate 
time  data  used  to  compare  the  original  signals  to  the  randomized 
signals.  If  the  randomized  signal  diverges  from  the  original 
signal,  this  signifies  that  the  time  series  data  is  nonlinear. 
[0037]  According  to  the  alternate  seventh  step  of  the  flow 
chart  70b,  the  results  of  the  nonlinear  analysis  of  the  time 
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series  data  are  compared  to  the  baseline  values  obtained  in 
second  step  45a  or  45b  to  approximate  the  damage  level  if  any  to 
the  structure.  The  level  of  damage  is  an  approximation  based  on 
the  differences  between  the  time  series  data  and  the  baseline 
values.  The  magnitude  of  the  difference  is  directly  related  to 
acceptable  tolerances,  which  are  based  on  the  application  of  the 
structure  itself.  A  static  structure  may  allow  for  more  of  a 
difference  between  the  data  and  baseline  values  than  a  dynamic 
structure . 

[0038]  According  to  the  alternate  eighth  step  of  the  flow 
chart  75b  a  decision  must  be  made  as  to  whether  to  generate  an 
alert  because  the  structure  is  seriously  damaged.  If  there  is  a 
serious  level  of  damage  as  indicated  by  a  large  difference 
beyond  the  acceptable  baselines,  then  the  graphical  user 
interface  will  create  a  visual  alert  80b  for  a  user.  If  there 
is  not  a  serious  level  of  damage,  the  next  set  of  time  series 
data  of  vibration  signals  is  analyzed  in  the  final  step  85b. 
[0039]  The  signal  processing  fault  detection  system  10  of  the 
present  invention  has  been  tested  on  a  helicopter  flexbeam  as 
illustrated  in  FIG.  3.  The  flexbeam  is  a  solid  metal  support 
member  of  a  helicopter  rotor  system.  It  has  a  bifurcated  end 
with  an  upper  appendage  and  a  lower  appendage.  As  part  of  the 
test  the  otherwise  intact  flexbeam  was  deliberately  damaged  in 
the  upper  appendage  to  introduce  a  fault  in  the  structure. 
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Strain  sensors  were  attached  to  both  appendages  while  the 
flexbeam  was  subjected  to  vibration.  The  vibrating  input  signal 
was  50kHz.  The  output  signals  from  the  strain  sensors  attached 
to  both  the  intact  appendage  and  the  damaged  appendage  are 
illustrated  in  the  time  series  graph  in  FIG.  4.  Within  the  time 
series  graph  shown  in  FIG.  4,  there  exist  time  delays  and 
amplitude  differences  that  are  qualitatively  observable. 

10040]  Mutual  information  analysis  is  performed  on  the  sensor 
output  signals  (i.e.,  time  series  data)  from  the  intact 
appendage  to  serve  as  the  baseline  values  for  a  "healthy" 
appendage,  corresponding  to  the  second  step  45b  of  the  method  in 
FIG.  2.  The  result  of  the  mutual  information  analysis  is 
illustrated  in  the  graph  in  FIG.  5.  Mutual  information  analysis 
is  then  performed  on  the  sensor  output  signals  from  the  damaged 
appendage,  the  results  of  which  are  illustrated  in  FIG.  6.  This 
corresponds  to  the  fourth  step  55  of  the  method  in  FIG.  2.  Both 
plots  of  mutual  information  indicate  linear  behavior  from  the 
upper  appendage  and  the  lower  appendage.  Linear  behavior,  as 
seen  in  the  graph  in  FIG.  5  is  verified  by  the  phase  randomized 
or  surrogate  data  not  diverging  appreciably  from  the  solid  line, 
which  represents  the  mutual  information  for  the  original  (i.e. 
non  randomized  phase)  waveform.  For  a  nonlinear  sensor  signal, 
the  surrogate  data  would  diverge  from  the  mutual  information 
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plot  for  the  reason  that  nonlinear  attributes  of  a  time  series 
are  contained  in  the  phase  portion  of  the  signal. 

[0041]  According  to  the  fifth  step  50  of  the  flow  chart  of 
the  method  in  FIG.  2,  the  data  is  determined  to  have  linear 
characteristics.  At  step  six  of  the  method  65a,  the  wavelet 
transmittance  function  { WTF )  ,  Equation  5,  is  selected  to  perform 
a  linear  analysis  on  the  data  obtained  from  both  sensors  on  both 
appendages . 


WTFl2M 


Km 

Km 


(5) 


WTFXt2[0,i)  is  the  ratio  of  sensor  1  average  Wavelet  transform  to 
sensor  2  average  Wavelet  transform.  The  Wavelet  transform  is  a 
2-D  function  of  frequency  o>  and  time  t  .  Wv(a),t)  is  the  average 

wavelet  transform  value  for  each  sensor  1,  and  2  respectively, 
which  transformed  into  the  frequency  arid  time  domains  using  a 
complex  wavelet  as  described  in  Equation  2.  A  graphical 
representation  of  the  wavelet  transmittance  function  of  the 
healthy  and  damaged  appendages  is  illustrated  in  FIG.  7.  The 
wavelet  transmittance  function  provides  a  method  to  quantify  the 
amount  of  damage  within  the  flexbeam.  This  corresponds  to  the 
seventh  step  70a  in  the  method  30  of  FIG.  2.  For  example,  an 
amplitude  level  of  1.2  may  indicate  an  excessive  amount  of 
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energy  in  any  appendage  of  the  fl exbeam  if  the  baseline  level  is 
set  at  an  amplitude  level  of  0.2. 

[0042]  In  another  test  of  the  fault  detection  system  10, 
piezoelectric  sensors  were  placed  on  two  aluminum  cantilever 
beams  as  illustrated  in  FIG.  8.  The  beam  dimensions  are 
l"xl^"x30" .  A  crack  was  deliberately  created  in  one  of  the 
cantilever  beams.  Both  the  damaged  and  the  undamaged  cantilever 
beams  were  clamped  to  a  test  fixture.  The  sensors  were  arranged 
on  each  beam  so  that  as  the  wave  front  propagates  across  the 
beam  the  crack  in  the  damaged  beam  is  between  the  sensors  as 
illustrated  in  FIG.  8.  One  sensor  measures  the  strain  created 
from  the  actuator  at  the  fixed  end.  The  second  sensor  measures 
the  strain  beyond  the  crack  on  the  damaged  cantilever.  The 
undamaged  beam  has  the  sensors  in  the  same  locations  so  as  to 
have  the  same  modal  effects.  Sensor  time  series  data  was 
obtained  for  both  cantilever  beams. 

[0043]  Mutual  information  analysis  was  performed  with  and 
without  random  phase  on  the  time  series  data  obtained  from  the 
sensors  attached  to  the  healthy  beam  to  obtain  baseline  values 
corresponding  to  the  second  step  45  in  the  method  30.  Mutual 
information  analysis  was  then  performed  with  and  without  random 
phase  on  the  time  series  data  obtained  from  the  sensors  attached 
to  the  damaged  beam.  The  results  of  the  mutual  information 
analysis  performed  on  the  time  series  data  corresponding  to  both 
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cantilever  beams  are  illustrated  in  the  graphs  in  PIG.  9,  which 
highlight  a  separation  or  divergence  in  the  plot  for  the  damaged 
beam.  Divergence  is  only  seen  in  the  time  series  data  obtained 
from  the  cracked  beam  because  the  crack  creates  a  nonlinear 
effect  in  the  data  and  thus  is  present  in  the  phase  randomized 
signal . 

[0044]  The  divergence  indicates  that  the  time  series  data  is 
nonlinear,  which  corresponds  to  the  fifth  step  60  of  the  method 
30.  Therefore,  a  nonlinear  analysis  is  performed  on  the  time 
series  data  according  to  the  sixth  step  65b  of  the  method  30 
using  mutual  information  analysis  combined  with  surrogate  data 
analysis.  The  results  of  this  combined  analysis  are  illustrated 
in  the  graph  in  FIG.  9.  A  determination  of  damage  level 
according  to  step  70b  of  the  method  30  is  accomplished  by 
measuring  the  divergence  in  the  two  mutual  information  plots. 
[0045]  The  advantage  of  this  invention  is  the  ability  to 
detect  nonlinear  behavior  in  complex  structures.  The  invention 
is  a  timesaving  analysis  tool  for  test  engineers,  which  could  be 
incorporated  in  a  function  generator.  Moreover,  this  tool  is  an 
aid  to  determining  the  optimal  signal  processing  technique 
required  to  analyze  sensor  data. 

[0046]  While  it  is  apparent  that  the  illustrative  embodiments 
of  the  invention  disclosed  herein  fulfill  the  objectives  of  the 
present  invention,  it  is  appreciated  that  numerous  modifications 
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and  other  embodiments  may  be  devised  by  those  skilled  in  the 
art.  Additionally,  feature (s)  and/or  element (s)  from  any 
embodiment  may  be  used  singly  or  in  combination  with  other 
embodiment (s) .  Therefore,  it  will  be  understood  that  the 
appended  claims  are  intended  to  cover  all  such  modifications  and 
embodiments,  which  would  come  within  the  spirit  and  scope  of  the 
present  invention. 
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ABSTRACT  OF  THE  DISCLOSURE 

A  fault  detection  system  designed  to  evaluate  the 
structural  integrity  of  a  material  employs  an  array  of  sensors 
disposed  over  the  material  being  evaluated.  The  sensors  detect 
vibrations  in  the  material  and  the  sensor  signals  are  fed  to  a 
data  processor.  The  processor  employs  a  method  to  analyze  the 
linear  and  nonlinear  characteristics  of  the  sensor  signals  and 
then  determines  whether  to  proceed  with  a  linear  signal 
processing  analysis  or  a  nonlinear  signal  processing  analysis  of 
the  sensor  signals.  Once  the  analysis  is  completed,  the  results 
are  compared  to  baseline  results  to  determine  what  if  any 
divergence  exists  between  the  results  and  the  baseline  results. 

A  significant  divergence  indicates  a  potential  material  failure. 
The  fault  detection  system  will  indicate  such  a  potential 
failure  through  a  visual  alarm  on  a  graphical  user  interface. 
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